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INTRODUCTION 


With the launch of the first Earth Resources Technology Satellite 
(ERTS)* on July 23, 1972, a new era of direct application of space 
technology to human needs began. Since that time, the usefulness of the 
mul tispectral image data obtained from the Landsat series of satellites has 
been demonstrated for aiding natural resource management. The ability to 
obtain synoptic measurements of systems spread over wide geographic areas 
offers an advantage over many traditional data taking approaches. Other 
advantages include easy access to data over 'inaccessible regions and 
economical frequent updates. Recent studies and programs such as Sigma 
Squared^iy, Large Area Crop Inventory Experiment ^ 2 ) • and others 
indicate remotely sensed data will be useful for agricultural forecasting. 


As experiments and studies such as Optimum Repeat Cycle Analysis 


(7) 


are identifying and quantifying the information obtainable via remote 
sensing, overall studies are nei^eded to address an operational system. 
Operational system studies must consider the data volumes, the computational 
times required, the logistics of gathering, processing, and disseminating 
the data, and the cost of impWmentation. In FY-77, the Marshall Space Flight 
Center Data Management Program activities centered around the analyses of the 
far term (1985 and beyond) Office of Application objective of applying space 
technology to an operational Global Crop Production Forecasting System. The 
Global Crop Production Forecasting Trade Study performed as part of this 
program identified major areas of disparity between projected technology and 
requirements. Out of this study, the concept of obtaining a repeated 
number of observations of sample regions via satellite was developed. The 
concept of reducing the data volume at the earliest point in the system was 
aljso advocated. From the previous studies, it is evident that remotely 
sensed multispectral image data will be a valuable and essential part of 
an operational agricultural system. However, there is no clearly decisive 
set of requirements currently established. Other studies such as Temporal 
TnWstigation for Mission Evaluatidtn ^gy sfe directed toward further 

defining the requirements as a function of the agricultural science. This 
study is directed toward establishing the data system relationship between 


*The ERTS was later renamed Landsat. 


the requirements and methods of implementation. Several questions 
concerning the economic trade-offs for an operational system were singled 
out for further definition in this study. These questions are in two major 
areas. One is the number of satellites required to reliably obtain 
cloud free images at precise times during the growing season. The other is 
tb better quantify the volume of data to be processed as a first step in 
optimizing a processing system. This report describes the investigation of 
these two factors. 


STATEMENT OF THE PROBLEM 

The goal of this study is to identify cost effective approaches to 
meeting the Global Crop Production Forecasting Objective. This implicitly 
favors approaches with a minimum number of satellites and a minimum amount 
of data processing requirements. 

Satellites are major fixed cost elements in an operational system. As 
such, it is desirable to design a system that will use the minimum number of 
satellites and still perform satisfactorily, it is also desirable that the 
amount of data processing required be minimized. Because of the temporal 
nature of the information in a cropland image, the positioning of a 
satellite as a function of time is a principal requirement. The main driver 
that must be satis ifed, which influences the minimum number of satellites 
to render satisfactory performance, is a function of this time dependent 
positional requirement and the number of picture elements needed to cover 
the world's agricultural areas. This number is a function of the spatial 
resolution required to obtain the desired accuracy. The image data acquired 
tising satellites, as depicted in Figure 1, is processed and combined with ^ 
other information called collateral data; 1) to measure or inventory the 
amount of land in production for particular crops, 2) to determine plant 
vigor as an indication of growth stage and potential yield, and 3) to 
assess the extent of stress from either environmental or induced episodes 
ajs it affects yield. Each of these uses of remotely sensed image information 
ijmposes constraints on the timing of df-ta acquisition. Of the three, the 
iinventory function is the best defined. Based on current United States crop 
production estimating practice, it is also the function most amenable to near 
term improvements from satellite sensed image data* 

2 




F I GURE-1 Tf5E-USE-0FL^ATELL l-TES-FOR AGRICULTURAL FORECASTING 



, This study is limited to accomplishing this inventory function only. 

It concentrates on determining the dependence of the success of obtaining 
the regulred Information on orbital altitude, period, position, and sensor 
swath. 

Other factors that indirectly impact the number of satellites needed and 
the data volumes encountered were considered. They form the baseline conditions 
for this study and are summarized in Table 1. The principal rationale for 
t^is baseline is that It is a reasonable set of conditions based on present 
technology and related studies. Any system expected to be operational in 
the 1985 time frame will be based on such current technology. 

‘A. 

: \\ 

CONCLUSIONS 

Results of this study support the following conclusions: 

1. Use two satellites of the Landsat-D design to accomplish the 
' agricultural inventory objective. 

I.'- The nominal Landsat-D will perform satisfactorily. 

3 . Accept scenes with up to S0% cloud cover for preprocessing. 

A. Use a regionally compensated oversampling strategy. 

5 . Extract samples Immediately after preprocessing. 

1" ' ' 

Additional results from this study indicate sufficient potential 
benefits to justify future study of the following alternatives: 

1. Use a floating sample approach with a reduction in acceptable 
scene cloud cover to as little as 30%. 

■'|2. Use dedicated satellites with on-bdafd sample extraction. 

3 . Increase the sensor swath width. 


TABLE! 

TRADE STUDY BASELINE CONDITIONS 





BASELINE 

VALUE 

RATIONALE 

Spectral Bands 

5 Visible + 1 IR 

Existing Thematic Mapper Design 

Spectral Resolution 

128 Levels 

Existing Thematic Mapper Design 

Spatial Resolution 

30 Meter Visible 
120: Meter IR 

Existing Thematic Mapper Design 

Swath Width 

100 NM 

Existing Thematic Mapper Des i gn 

Crops 

Wheat, Corn, Soybean, Rice - — 

Trade Study 1 Results 

Countries 

Regions 

22 ' 
Major Ecological 

Contribution to World Crop Value 
ECO Systems 

Sun Synchronous 
Increased Swath 

Yes - ^ 

Minimum Change 

Analysis 

1985 Operational Need 

Low Altitude 

Grazing Angle 

1985 Operational Need , Landsat 1 ,2 
Experience; BRP 

Window Length ^ 

Discrimination 

■ ■ ' j ■ ■ 

ECO Systems (Ref. 8) 

Confidence Limits 

1/20 Year 

ECO Systems (Ref. 8) 

Sample Segment Size 

1 KM 

Discrimination 6 Confidence 
ECO System (Ref. 8) ^ 

No. of Samples 

1000 Base 

ECO Systems (Ref. 8) 

Use of Samples 
Cloud Condition 

Yes 

Trade Study 1 (Ref. lO) 
Allied Cloud Model 



STUDY APPROACH 


j The approach used for this investigation was to model the satellites, 
their movement and observation capability, world cloud conditions, and major 
growing regions for wheat, corn, soybeans, and rice using the general simuia- 
jion capability of MSFC's Data System Dynamic Simulator (DSDS)^j^j. Twenty- 
two countries were chosen based on the criteria that each contributed two 
percent or more to the world harvest of one of the selected crops. The larger 
countries were divided into geographic regions corresponding to statistical 
reporting districts. The point target capability of DSDS was used to provide 
a statistically sound measure of variations due to geographic location and 
cloud conditions. A minimum number of samples in a simulation region was set 
at 30 for a geog raphi e region containing on iy one c rop. The number of samp 1 es 
was adjusted to a maximum of 60 when all four crops were grown in a region. 

I This approach resulted in 1553 samples for the world (see Table 2) and 
(j»ermitted the collection of statistics on individual samples and regions. This 
permitted a measure of the effects of opportunity variations due to overlap in 
coverage, length of time during which observations could be obtained, and the 
effects of localized cloud conditions. The number of samples drawn for an 
operational system varies for each region according to local accuracy con- 
siderations and can be scaled from the simulation data. Table 2 lists each 
country, its regions and crops, the number of samples used in the simulation 
model and the number of samples estimated to be needed for an operational 
system. 

I - Alternatives to optimizing the number of satellites were considered. A 
working Blue Ribbon Panel", made up of technical people working in related 
disciplines, considered the technical feasibility, the practicality, the cost 
impact and any additional constraints likely to be encountered for each 
alternative. Those alternatives appearing most likely to be successful were 
Conceptualized as candidate systems. The Data Systems Dynamic Simulator was 
usied to measure parametric variations in the candidate systems. 

*This panel consisted of General Electric and NASA personnel and convened for 
the purpose of this Trade Study. The panel met in Beltsville, MD, Nov. 9, 1977 
Individual panel menbers subsequently provided consultation for this study. 


TABLE 2. COUNTRIES, REGIONS, CROPS, AND NUMBER OF SAMPLES USED FOR STUDY 


Country & Region 


Argentina 
Australia 
Bangladesh 
Brazii North 
Brazil South 
Canada 
China North 
China Central 
China South 
Egypt 
France 

India Punjab 

India Ganges ' 

India Central 

India Bilaspur 

India Coastal 

Indonesia 

Italy 

Japan 

Mexico 

Pakistan 

Romania 

South Africa 

Phi 1 ippines 

Thai land 

Turkey 

USA - Region A 
USA - Region B 
USA - Region C 
USA - Region D 
USSR Latvia 
USSR Ukraine 
USSR Transvolga 
USSR Voiga-Ural 
USSR Siberia 
Yugoslavia 

TOTAL 


Crops 


W,C 

W 

R 

C 

C,S,R 

W 

W,C,S,R 

W,C,S,R 

W,C,S,R 

C 

W,C 

W,C 

W,C,R 

W,C,R 

W,C,R 

R 

R ) 

W,C,R 

R 

C 

W 

W,C 

C 

R 

R 

W 

W,C,S 

W.C 

W,C,S,R 

W,C,S,R 

W,C 

W.C 

W,C 

W,C 

W.C 

W,c 


Number of 
Simulation Sam 



Estimate of H of Oper- 
ational Sample Segments 


1500 

1000 

1000 

1000 

1767 

1000 

2000 

2000 

2000 

1000 

1500 

1500 

1767 

1767 

1767 

1000 

1000 

1767 

1000 

1000 

1000 

1500 

1000 

1000 - 
1000 
1000 
1500 
2000 

! 2000 

2000 
1500 
1500 

1 1500 

1500 
1500 
1500 


51355 


Wheat , C * Corn, S » Soybeans, R 
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ALTERNATIVES INVESTIGATED 


Some alternatives for optimizing the number of satellites were investigated 
and a brief disposition of each is included in Appendix A. Four candidate 
systems were chosen for further exploration using DSOS. They all have the 
property of using the currently designed Landsat-D vehitle, in some cases 
with modification to the sensor system or operational procedures. The four 
systems with their characterizing features are illustrated in Figure 2. Except 
for System 2, the orbital altitudes are sufficiently close to the nominal 
704 KM design altitude of the Thematic Mapper scheduled for use on Landsat-D • 

I 

that a sensor and satellite redesign will not be required. As King^l2j 
described, an operational system could use different altitudes and satellite ^ 
pals it ions at different times according to needs with a small impact on the 
consumable budget. The four candidate systems compTise combinations of 
satellites in four different sun'-synchronous orbits. Each of these orbits were o 
investigated using OSDS. 


704 KM ORBIT - NOMINAL LANDSAT-P 

I A sun-synchronous orbit at an altitude of 704.052 KM and an inclination - 
of 98.204 degrees provides total earth coverage every 16 days. This is the 
design orbit of the Landsat-D with the Thematic Mapper covering l85 KM swath 
on each orbit. Every 233 orbits, thie nadir repeats itself. 


It has the advantage of not requiring special programming of its position 
as well as providing complete global coverage. It is inefficient from an 
agricultural viewpoint because it is often over uninteresting regions such as 
oceans or deserts. The time between repeat coverage of 16 days is unacceptable 
af certain times either because of the temporal nature of the information or 
blecause of delay in acquiring information unobtained due to atmospheric 
conditions. To overcome these difficulties, additional satellites may be 
used. One outcome of the Study was a measure of the effectiveness of 
additional satell ites at this Orbit versus other orbits for obtaining the 
desired samples. 


CANDIDATE SYSTEM I 



CANDIDATE SYSTEM 3 



CANDIDATE SYSTEM 2 


TWO SATELLITES 1485 KM EACH 20 DAY REPEAT — 


185 KM SWATH - 704 KM alt. 


NO. OF TOTAL WORLD 

satellites COVERAGE cycle 

1 ' 16 DAYS 

2 8 DAYS 

3 5 1/3' DAYS 


! LANDSAT-D S ATELL ITE 


CYCLE INCREASED SWATH 



TOTAL WOR LD COVERAGE E VERY _ 
10 DAYS (OR LESS) AS 
FUNCTION OF SWATH. 


MOO I F I CD LANOSAT-D SATELL 1 TE 


CANDIDATE SYSTEM 4 


• ' ^ ' i I * 

ONE SATEUJJ'E_JL_70k^^^^ 
DAY REPEAT)' 

ONE SATELLITE AT 743 KM 
(9 daV repeat) 


LANDSAT-D SATELLITE 



ONE SAfELLITE 704 km 
(16 DAY) 

ONE SATELLITE AT 725 KM 
(2 DAY) 

(PROGRAMMABLE POINTED) 

LANDSAT-D SATELLITE 
GROUND SCHEDULING SYSTEM 


FIGURE 2. FOUR CANDIDATE AGRICULTURAL DATA COLLECTION SYSTEMS 


7A3 km orbit 


■ III,, .1 ■■i-i ^ . . I ■ 

An increase of 39 KM in altitude provides an orbit that repeats after 130 
orbits. This is insufficient to provide total earth coverage with I85 KM 
swath. It covers only 2A,050 KM or about 60 percent of the area at the equator. 
It does offer the advantage of 9-day repetition, which when combined with other 
orbits has the potential to provide optimum coverage of areas of agricultural 
interest. Additional capability of an orbit that repeats frequently can be 
obtained using the stabilization system in the Landsat-D which can be biased 
to provide some off nadir pointing. One swath width is the usually accepted 
limit. Because of these possibilities, this orbit of 7^2.57 KM and a sun- 
synchronous inclination of 98.367 degrees was investigated. 

725 KM ORBIT 


1 A potentially useful orbit between the 70k KM 16-day repeat and the 
7^3 KM 9“day repeat occurs at 72i».35 KM. It repeats itself every 29 orbits or 
two days. While it only covers about i/8 of the world at the equator, it has 
the advantage of a short repeat cycle. When combined with satellites at other 
oi'bits, and use is made of one swath pointing as well as position adjustments, 
areas of specific Interest could be frequently observed. Such a system assumes 
the dedication of the satellite to the needs of the agricultural user so it 
can be pointed or positioned without concern for loss of data to other users. 
The sun-synchronous Inclination for this orbit is 98.292 degrees. 



1A85 KM ORBIT 

! One additional orbit was investigated at 11(84.65 KM, with a sun-synchronous 
Inclination of 101.874 degrees. This altitude, being about twice the normal 
Landsat-D orbit, would require a redesign of the Thematic Mapper^ 


^ While repeating at a longer period, every 20 days, 1 1 might pertnl t greater 
flexibility of pointing without encountering unacceptable radiometric and 
geometric difficulties due to the grazing angle with the earth. This grazing 
angle Is discussed In Appendix B. s~, 


INVESTIGATION PROCESS 


The Data System Dynamic Simulator (DSDS) at MSFC was used to determine 
the foTlowIng: 

o Orbital position of the sateilltes as a function of time and orbital 
parameters. - ','7' 

o The sensor field of view as a function of satellite position and 
swath width. 


o The state of the samples within the field of- view as a function 
of the cloud model. 

o The statistics on the target acquistions as a funetlon of the 
above, the preprocessing acceptance criteria and the processing 
requirements. 


The DSDS combined the dynamics of satellite position, crop models, 
clloud models, and processing requirements. A Monte-Carlo method was used in 
conjunction y/lth a cloud model to Inject the realism of cloud cover Into 
the simulation. 


The simulation was segmented as Illustrated in Figure 3« For economy 
of simulation time, the results of each successive simulation Segment were 
saved and used for later parametric variations. For example, for a one year 
simulation of satel 1 i te positions repeating , every 16 days, the mission 
ephemerls was generated for 16 days and reused. For different sensor swath 
widths for the same satellite positions, the same mission model data was 
lised. The results of the Mission Model, sensor swath, and target described the 
data available from one satellite. Any combination of satellites desired were 
combined in the multivehicle crop model to constitute one of the candidate 
systems. Thus, such subtle effects a® the difference In insertion tim and 
position offset for a 3 Landsat vehlc^^ system were fncluded. The effect of 





o Combine mu] tip 
sate! 1 i tes 




cloud conditions were determined by a comparison of a random number against 
the predesignated cloud condition breakpoints for the cloud region correspond- 
ing to each sample. The acceptable cloud conditions for preprocessing was 
the first criterion to be met for all the samples in a scene 90 nautical miles 
aljOng track and one swath width wide. A subsequent random number was com- 
pajred with the scene cloud conditions to determine if each sample was clear or 
cloudy. A more detailed discussion of the DSDS models and the available data 
Is presented in Appendix C. 


PERFORMANCE OF CANDIDATE SYSTEM 1 WITH 1. 2, OR 3 SATELLITES 


j A major portion of this study consisted of establishing the capabi 1 Ity of 

■' ‘ ^ ^ 

obtaining cloud free observations of the desired sample segments with nominal 
Landsat-D satellites. 


The bar chart In Figure gives a rapid comparison of the results achieved 


with 1, 2, and 3 

satellites. Each bar In Figure k represents the mean percent 

of samples acquired at the 50% confidence level for each 

of the 36 regions. The 

last bar for each of the 3 cases 

“is the world average. - 

The data Is taken from 

columns 1, 5> and 

9 of Table 3. 

Additional Information 

from Figure 4 Is summar- 

ized below. 




i 

Mean Percent of 

i ~ 

Percent of Regions 

Minimum Percent 

Number of 

Desired Obs. 

Obtaining 98% 

Achieved for 

Satell i tes 

Achieved 

of Needed Samples 

Any Region 

;t:' 

87.1 

■ -.li.e 

,1 . 

72.3 

: 2 

97.0 


31.6 

3 

99.2 

88.8 

96.3 


The mean and standard deviation for the percent of samples for which 
loud free observations were missed are presented in Table The two satellite 
data Is based on 20 simulation runs and the three satellite data Is based on 
7 runs. The 50% and 95% confidence limits for obtaining samples were deter” 
mined for this data. The 95% confidence level corresponds to cloud free 


*As a vehicle of distinction, the term sample segment will be used when referring 
to the image data 1 KM x I KM in the modeled dperatlonal system and the Word 
sample will be used when referring to the point repreSen^^^ the simulation:. 
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PERCENT OF SAMPLES ACQUIRED 



1 SATELLITE 2 SATELLITES 3 SATELLITES 

FIGURE k. COMPARISON OF 1, 2 , AND 3 LANDSAT-D SATELLITES 


Note: Each bar represents a region. 


ORIGINAL PAO* IK 
OP POOR QUALITY 





TABLE 3. PERCENT OF SAMPLES ACQUIRED FOR EACH TEST CASE 


COLUMN NUNSCR 


NUMSCR Of VEHICLES 


VEHICLE ALTITUOE 

SWATH WIDTH (KM) 

MAXIMUM CLOUD 
COVEA ACCEPTED (t) 


NO. Of AUNS AVEAA6E0 

WOALO AVEAACE 

AAC^NTINA 
AUSTAALIA 
SANCLAOESH 
BAAZIL - NOATH 
BAAZIL - SOUTH 
CANADA 

CHINA - NOATH 
CHINA - CEHTAAL 
CHINA - SOUTH 
EGYPT 
fAANCE 

INDIA - PUNJAS 

INDIA - GANCES 

INDIA - CEHTAAL 

INDIA - BILASPUA 

INDIA - COASTAL 

INDONESIA 

ITALY 

JAPAN 

MEXICO 

PAKISTAN 

KOMANIA 

S. AfAICA 

PHILIPPINES 

THAILAND 

TUAREY 

US - A 

US - B 

US - C 

US - D 

USSR - LATVIA 
USSR - UKRAINE 
USSR - TRANS-VOLGA 
USSR • V0L6A-URAL 
USSR - SIBERIA 
YUGOSLAVIA 



















































































ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE MEAN AND STANDARD DEVIATION FOR THE PERCENT 
OF SAMPLES MISSED* 


COUNTRY/REGION 


EGYPT 

AUSTRALIA 

CANADA 

THAILAND 

USSR - LATVIA 

USSR- VOLTA-URAL 

USSR - UKRAINE 

USSR - TRANS-VOLGA 

USSR - SIBERIA 

BRAZIL - SOUTH 

BRAZIL - NORTH 

JAPAN 

ARGENTINA 

MEXICO 

S. AFRICA 

ROMANIA 

ITALY 

FRANCE 

TURKEY 

PAKISTAN 

YUGOSLAVIA 

PHILIPPINES 

INDONESIA 

USA - REGION B 

USA - REGION C 

USA - REGION A 

USA - REGION D 

CHINA - SOUTH 

CHINA - NORTH 

CHINA - CENTRAL 

BANGLADESH 

INDIA - PUNJAB 

INDIA - COASTAL 

INDIA - BILASPUR 

INDIA - GANGES 

INDIA - CENTRAL 

WORLD 


2 LANDSAT-D VEHICLES 


MEAN 

STD. DEV 

0 

0 

0.17 

0.74 

0.17 

0.74 

0.17 

0.74 

0.06 

0.25 

0.06 

0.25 

0.17 

0.41 

0.17 

0.54 

0.78 

1.30 

0.09 

0.42 

0.50 

1.22 

0.33 

1.02 

O.U 

0.66 

O.kk 

0.66 

0.67 

1.74 

1.17 

1.22 

1.26 

1.10 

1.50 

1.41 

1.50 

2.29 

1.67 

2.54 

2.00 

2.06 

2.25 

2.61 

3.4i» 

2.52 

1.13 

1.18 

3.88 

1.15 

i».65 

2.26 

5.15 

2.26 

1.25 

1.70 

2.64 

3.57 

8.42 

2.27 

5.66 

5.72 

0.78 

1.57 

2.16 

1.51 

5.19 

3.67 

8.07 

6.61 

8.30 

4.40 

3.01 

0.50 


3 LANDSAT-D VEHICLES 


MEAN STD. DEV. 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.32 

.84 

0 

0 

0 

0 

0 

0 

0 

• 0 

0 

0 

0 

0 

0 

0 

.13 

.35 

.16 

.42 

0 

0 

0 

0 

.32 

.54 

.64 

.87 

.64 

.87 

0 

0 

.40 

.50 

2.11 

.99 

3.50 

1.16 

0 

0 

.83 

.83 

1.90 

1.15 

.32 

.54 

0 

0 

.48 

.59 

3.69 

3.77 

1.14 

.81 

2.43 

1.97 

.82 

.30 


*Fai1ure to obtain 1 cloud free observation during the designated time 
for each sample segment constitutes a miss. 
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observation in 19 out of 20 years. The cloud free samples required for the 
36 regions are shown in Figures 5 and 6. Solid bars represent 95% confidence 
level and the crosshatched bars represent 50% confidence level. On Figure 5 
(two satellite case), it can be seen that only 5 of the 36 regions will 
furnish less than 90% of the sample segments. On Figure 6, it can be seen 
that with three satellites, none of the 36 regions will have less than 90% 
of the sample segment needed. These resul ts are based on only using scenes 
with 50 % cloud coyer or less. 


COMPARISON OF CANDIDATE SYSTEMS 


The candidate systems were compared based on their capability of obtain- 
ing 98 % of the designated Samples in each region. Each candidate system is 
described in Figure 2. 


The percent of samples acquired regionally for each test case is presented 

j - ■■ 

|n Table 3. Test cases 1 through 9 apply to candidate system 1; Test cases 

10, 11, and 12 apply to candidate system 2; and Test cases l^i and I 3 apply to 
candidate systems 3 and 4, respectively. A summary of this data is shown in 

Figure 7 . 


Of the four two-satellite candidate systems, the 70i» KM orbit or nominal 
iiandsat-D performed the best. The combination of one vehicle at 704 KM and one 
at 7^3 KM, system 3» performed nearly as well and has the advantage of providing 
additional viewing opportunities when repeat coverage is required. System 2, 
vjfith two vehicles at 1485 KM, performed approximately the same as System 3. 

A disadvantage of System 2 is that a redesign of the sensor would be required 

1 

because of the substantial increase in altitude from the nominal Landsat-D. 
System 4 performed poorer than any of the other candidates. This is because 
the two-day repeat satel 1 1 te at 725 KM on has an opportunity to view 20.3% 
of the needed samples. 
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EFFECT OF INCREASED SWATH WIDTH 


Comparison runs for a 20 percent increase in swath width were made for 
candidate systems 1 and 2. The percent of samples acquired, the increase in 
samples acquired due to increased swath width, and the reduction in the number 
of samples missed for each case are tabulated below: 


TABLE 5. EFFECT OF 20 % INCREASED SWATH WIDTH 


Number of 
Satel 1 i tes 

Altitude 

KM 

Samples 
- Perc 

Acqui red 
.ent 

Samples 
Acqu ired 
% Increase 

Samples 
Missed 
% Reduction 

185 KM Swath 

222 KM Swath 

1 

704 

87.1 

92.3 

6.0 

40.3 

2 

704 

97.0 

97.9 

0.9 

30.0 

1 

1485 

■ 83.1 

88.9 

7.0 

34.3 

2 

1485 

96.5 

97.8 

1.3 

37.1 


I As can be seen, a 20 percent increase in swath width causes a maximum of 
7.0 percent increase in the samples acquired. However, there is a reduction in 
the percent of samples missed of 30*0 to AO. 3 percent. Thus, a small percentage 
increase in swath width yields approximately a 2 to I percentage reduction 
In the number of samples missed. 

The effect of increased swath width is shown in Figure 8 for one Landsat 
Satellite and in Figure 9 for two satel 1 ites. In both figures, the percent of 

samples acquired in each of the 36 regions is shown for four cases. Starting 

on the left side of the figure, the four cases are: — 

! . . ' 

o j Nominal Landsat-D 

o '20 % cloud cover accepted - 

' ^ o 90^ cloud cover accepted 

o increased swath width and clouds 

F 8 and 9 give a quick appreciation for the benefits of increased 

Cloud cover acceptance and swath width. Two satellites with increased swath and 
cloud cover acceptance perform si ightly better than the nominal three satellite 
case. 
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EFFECT OF INCREASED CLOUD COVER ACCEPTANCE 


The effect of accepting scenes with up to 90% cloud cover upon the 
abi 1 i ty to obtain needed samples was investigated for one and two Landsat-D 

'V' ’’f ■; " 

satellites. The results from the four cases compared are shown beTow: 


Number of 
Satell i tes 

Maximum Cloud 
Cover Accepted 
{Percent) 

Samples Acquired 
(Percent) 

Samples 
Acqu i red 
% Increase 

Samples 
Missed 
% Reduction 

1 

0 : 

50 

87.11 


^ ■■■' : 1 ^ 

; - ; 1 

1 

90 

92 . 11 » 

5.8 

39.0 

2 

50 

96.99 



2 

90 

98.8I1 


61.5 

1 Increasing the acceptable 

cloud cover at the 

preprocessor 

had a significant 

e'ffect on 

reducing the percent 

of samples missed 

for the two satellite cases. 


f 


0,n a global basis, it reduced the percent of samples missed by 61.5%. For the 
ive regions in Figure 5 which did not achieve 90% of the needed samples at the 
9|5% confidence level, the mean percent of samples acquired increased from 92.86% 
to 96.7^%. Thus, through regionally selected use of scehes with up to 90% 
cloud cover and a small -percentage of dve r samp I log (up to 10%)', the 98% of 
samples needed for each region would be obtainable. * ' 


The effect of increased cloud cover on the processing load is shown below; 



Maximum 


NUMBER OF 

SCENES TO 

BE PROCESSED 

Cloud 

Category 

Accepted 

Cloud Cover 
Accepted 
(Percent) 

Samples 
Acqu i red 
(Percent) 

Total 
For Year 

Peak 

Day 

Average for 
Peak Period 

2 

1 1 

' 3P 

1 1 

9f».09 

1»790 

ki 

36 (26) 

3 

. 50 

96.99 

6383 

59 

k8 (3k) 

k 

90 

38. 8k 

93.17 85 70 (50) 

5 Day Work Week-^ ^ 


(7 Day Work Week) 


2k 


The average number of scenes to be processed per day during the peak 
growing season from May to September is based on processing only one obser- 
vation per window. This was chosen to measure the sensitivity of obtaining 
samples versus the number of sate1tites\ The actual number of observations 
requi red _is expected- to be greater^ but the exac% number has not been 
determined. 


From Figure 5, it can be seen that in order to obtain 98% of needed 
samples, scenes with up to 90% cloud cover would be required in 25 of the 
36 regions if minimum oversampling is used. If oversampling is used to 
obtain all additional samples in regions where at least 92% of the samples 
are obtained, then scenes with up to 90% cloud cover would only be required 
in 9 of the 36 regions. Thus, depending upon the amount of oversampling 
used, the average processing load to obtain one observation per window 
during the peak period will fall somewhere between and 70 scenes .per day. 


EFFECT OF WINDOW LENGTH 


I The primary- factor effecting a given system's capabil i ty of obtaining, 
needed sample segments is the observation period or window length. In 
this study, window lengths were set to maximize the probability of discrlmi-l 
nating each crop from the.,other major crops and confusion crops in a given 
phenological regfon. Window lengths varied from 18 to 76 days (11/8 to 
k 3A repeat cycles). 


The samples in the United States were divided into four sets based on 
window length. The mean and standard deviation of the samples for which 
clear observations were missed are listed below as a function of window 
length. This data is based on a 20 year simulatiqn for a system with twQ 


Landsat-D satel 1 i tes> 


Set 

Range of 
Window 
liength 
1 (Days) 

Average 
Wi ndow 
Length 
1 (iiiays) 

Samples 

'Mean 

Missed (Pe 
Standard 

rcent) 

Deviation 

1 

1 18 

- 18.0 


4 

.14 

2 

28 

28.0 

4.75 

1 

.78 

3 

32-<i7 

41.7 

0.61 

0.70 


65-76 

68.0 

0.06 

0 

.26 




25 


I The 50 % and 35 % confidence limits for the samples acquired as a 
fLnct ion of window length are plotted in Figure 10. A window length of 
40 days is required to obtain 98% of the needed targets at a 35 % confidence 

I level. 

y 


The three Landsat-D satellite system was run with fixed window lengths 
of 15, 20, 25, and 30 days. The results plotted in Figure 11 are for the 
U.S. alone and for the average of all regions. Samples missed in the U.S. 
was below 2% for a 24-day window and below 2 % for the world for a 28~day 
window. 


RELATIVE OBSERVATION PERIOD 


In addition to window length, the overlap between adjacent swaths 
affects a system's capability of obtaining clear observations. For each 
rdgion. Table 6 lists the meaa window length, the mean latitude, overlap, 
relative observation period and percent of targets missed. The relative 
observation period is calculated as: Relative Observation Period = Window 

Length X (1 + Overlap). The percent of targets missed in each region is 
plotted in Figure 12 as a function of the relative observation period. 

A curve, drawn through the points, crosses the 1 % line with a relative 
observation period of 56 days. The large spread about the curve is due to 
differing cloud conditions for regions with the same relative observation 
period. 


MULTIPLE OBSERVATION PERFORMANeE 


To determine the capabll ity of obtaining more than one observation 
during a window, the number of observations accepted for processing was 
varied from 1 to 6 per sample. 


■ Figure 13 shows the capability of obtaining multiple clear observations 
from one satellite for each of the four orbits Investigated. This shows the 
advantage of 725 KM orbit If frequent coverage of a limited area Is required 
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TABLE 6. RELATIVE OBSERVATION PERIOD PER REGION 


I 



Country/ 

Region 

Targets 

Missed 

(Percent) 

Mean 
Lati tude 

Overlap 

■SSmilll 

Relative 
Observation 
Per 1 od 

(Deg rees) 

Mean 


Egypt 

0 

26.6 

0.203 


63.0 

' ^5.8 


Australia 

0.17 

■ 

-24.5 

0.182 


54.0 

63.8 


Canada 

0.17 

57.7 

1.002 


54. 0 

108.1 


Thai land 

0.17 

12.3 

O.lOi 


64.0 

70.5 


USSR - Latvia 

0.06 

58.0 

1 .030 

45.5 

92.4 


USSR - Volga-Ural 

0.06 

52.0 

0.747 

48.0 

83.9 

lUSSR - Ukraine 

0.17 

48.5 

0.624 

1 55.0 

89.3 

USSR - Trans-Volga 

0.17 

50.0 

0.674 

46.5 

77.8 

USSR - Siberia 

0. 78 


54.0 

0.830 

37.0 

67.7 

Brazil - South 

0.09 


-23.5 

0.173 

75.0 

88.0 

Brazil - North 

0.50 


-8.0 

0.086 

67.0 

72.8 

Japan 

0.33 


38.0 

0.365 


64.0 

87.4 

Argentina 

O.kk 


-32.6 

0.276 


62.0 

79.1 

Mexico 

OAk 


23.5 

0.173 

69.0 

80.9 

S. Africa 

0.67 


-28.3 

0.222 

i 62.0 

75.8 

Romania 

1.17 


46.0 

0.549 

45.0 

69.7 

Italy 

1.26 


42.0 

0.449 

45.0 

65.2 

France 

1.50 


47.0 

0.577 

45.0 

71.0 

Turkey 

1.50 


38.0 

0.365 

i 40.0 

54.6 

Pakistan 

1.67 1 

29.0 

0.230 

40.0 

49.2 

Yugoslavia 

2.00 

44.2 

0.501 

45.0 

67.5 

Phil ippines 

2.25 

11.0 
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f6r evaluating the extent or intensity of an episode. For total coverage, 
with one or two observations required during a window, the 704 KN orbit is 
clearly the best. 


I Similar plots are shown in Figure 14 for the 1, 2, and 3 satellites at 
7,b4 KM, and in Figure 15 for 1 and 2 satellites at 1485 KM. Figure l4 shows 
that a second clear observation can be obtained for 96% of the targets with 
three sateiiites deployed and for 87^ if two satellites are used. 

f.- ^ _ 

ORBIT INSERTION VARIATIONS 


Various insertion points were used for each of the candidate orbits to 
obtain equal time spacing between observations for the full coverage orbits 
(704 and 1485 KM) and to optimize coverage for the partial coverage orbits 
(725 and 743 km). As expected, different insertion points had hegiigibie 
effect on the capability of obtaining samples for the full coverage orbits. 


The insertion point had a significant effect on the percent of samples 
acquired by the 725 KM (2-day repeat) orbit. The percent of targets observed 
ranged from 18.2% to 20.4%, thus, the later insertion point viewed 12% more- . 
targets than the first. 

.. ■ * ■■ ■ ! ■ ■ ' . ^ 

EQUATORIAL CROSSING TIME 


1 

The current and planned Landsat orbits have a 9:30 AM equatorial crossing 
time. To determine if a later equatorial crossing would increase the proba- 
bility of achieving the desired samples, three simulation runs .were made with 
1:,00 PM equatorial crossing. The change in cloud conditions betWjeeii 9:30 AM 
and 1:00 PM did not have a significant effect. The inea-n percent desired 
saiDples acquired from the comparison runs is shown below: 


Equatori al Cross i ng Time 

■ i 

9:30 AM 
1:00 PH 


Samples Acquired (Percent of Desi red) 

mzk ■ 

99.28 









SUMMARY OF ORBITAL PARAMETER INVESTIGATION 


The investigation of the four candidtate systems included the 
following parametric variations: 

I 

o Orbital characteristics 

! 

o| Number of satellites 
o| Swath width 
o Cloud cover acceptance 
o Window length 
o Relative observation period 

I 

o Multiple observations 
o Orbital insertions 
o Equatorial crossing time 

The analysis to this point is based solely upon data acquisition without 
regard to data processing volume. 

For the purpose of obtaining areal measurement of agricultural land, 
candidate system 1, consisting of equally spaced satellites at 70A KM performed 
b|st. The quantitative relationship between the data acquisition and the 
various parameters were obtained. Data on the benefits of increased swath 
was determined and will be suitable for future comparisons of the costs and 

complexity of increasing the swath. Significant insight into the potential 

I 

benefits of altering the processing acceptance criteria was obtained. The 

dependence of performance upon window length and the number of observations 

V- ! 

required was explored. This data is valuable in extending the results 
of this study at the time- that the needs of the agricultural scientists for 

i ’ 

remotely sensed Image data are better defined. Data was also obtained ..that 
showed the results of this study are insensitive to orbital insertion and 
equatorial crossing time. 
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IMAGE DATA PROCESSING 


The data as acquired by the Thematic Mapper is a sequence of bits repre- 
sejnting discrete levels of energy received at different times in different 
pijedetermined energy bands. For efficient usage, the band and the pixel re- 
presented must be identified by an ordering of the data into a convenient format. 
Some pixels represent known energy levels and serve as reference points for 
caiibrating the data. Additionally, for agricultural usage, the pixels repre- 
sented must be identified to some reference system. Each of these functions, 
formatting, radiometric correction, and registration are considered preprocessing 
Tlj»e amount of computer resources required to perform these functions depends 
upon the complexity of the function, and the amount of data upon which the 
functions are performed. By employing a sample extraction philosophy, the 
aipount of data subjected to each succeeding process can be reduced. Part of 
this study was directed at obtaining realistic data loading at each functional 
point in the processing system. 


YEARLY LOADING FOR AGRICULTURAL USAGE 



The sizing and configuration of an operational processing system is 
dependent upon the time line of data acquisition. Two convenient measures of 
data volume are acquired scenes and the number of sample segments* A scene 
represents the 38,027,776 pixels as would comprise an area 185 KM squared at 
30 meters spatial resolution. Each pixel represents effectively 5.0625* bytes 
of data. A sample segment represents I 600 pixels in a desired area one KM square 
with kk% excess pixels to allow for registration error. A yearly plot repre- 
senting weekly scene and target acquisition is presented in Figure 16. This 
data was obtained for the condition of two satellites in a ,70^ KM orbit, 185 KM 
Swath, an acceptance criteria of 50^ or less cloud cover oh a scene, and the 
acceptance of only one observation per target per window. The targets plotted 
arc the simulation targets with a base number of 30 per single crop region. 

Each 30 targets represents a larger number of samples in an operational system. 

*Spcct ra 1 reso I u t i on is 256 levels or one byte/band for five bands at 30 meters 
spatial resolution plus one band at 120 meters resolution or (30/120)2 - .0625. 
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FIGURE 16. WEEKLY SCENE AND TARGET ACQUISITION FOR 1 YEAR 


The actual number will depend on the crop densities for the region. A figure 
of 1,000 samples Is a reasonable average and was used In determining processing 
loads. 

PROCESSING SYSTEM DATA VOLUMES 

The total number of pixels per year to be subjected to the differing 
processes was determined for the different candidate systems. This information 
Is shown In Figure 17> What Is readily apparent is the small percentage, less 
than 0.1^, of agricultural Information In the acquired data. This factor Indi- 
cates a need to extract th^ useful Information from the data at the earliest 
practical time In the sequence of agricultural data processing. 

An obvious way to reduce the data volume Is to reject Images of less than 
some established quality, as determined by the percentage of cloud cover In the 
scene. This Is the currently accepted operational approach. While It reduces 
the amount of unproductive activity of preprocessing data that will ultimately 
be rejected as unusable. It has the detrimental quality of also rejecting 
some highly desirable data. For the data shown in Figure 17. the acceptance 
criteria was established at 50%. The subsequent numbers represent data that 
could be expected In the 50% clear portion of the scenes accepted for prepro- 
cessing. Thus, the two columns of numbers; one headed by "Extract Agricultural 
Samples" which are all of the samples In the accepted scenes, and the other 
column headed "Edit Cloud Samples", which are only those that are In the clear 
part- of the accepted scenes. The latter column Is the yearly total In millions of 
pixels that can be expected for the different acquisition systems represented 
by each row In the figure. While the usefulness of multiple observation of the 
same sample segment Is acknowledged, this study placed a premium on the value 
of the first observation of a sample segment during each of the specified windows. 
Those observations are represented In the column headed "Process New Information". 
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♦MAXIMUM OF 50 PERCENT CLOUD COVER IN A SCENE ACCEPTED FOR PREPROCESSING. 


FIG URE 17 . A COMPARISON OF SATELLITE EFFECTS ON DATA VOLUME FOR PROCESSING 



COMPARISON OF CANDIDATE PROCESSING SYSTEMS 


A knowledge of the data volumes at each functional location is necessary 
tjo size the processing system. The required processing system will influence 
the choice of the data acquisition system. Nevertheless, the mission suita- 
bility of any given system must be measured in terms of its ability to 
acquire the required data. Because of regional variations discussed earlier, 
any single figure of merit such as those numbers in the column headed "Percent 
of Desired Samples Acquired" should be used cautiously. However, they do 
provide an additional dimension in addition to data volumes and the information 
to data ratio, by which the candidate acquisition system may be judged. A 
goal established early in the study was to acquire SB% of the desired samples. 
Only the three satellites system does that, although others come close. When 
the incremental cost of three versus two satel lites is taken in perspective, 
there is ample room to investigate other remedies. Oversampling is one 
serious consideration. 


OVERSAMPLING TO COMPENSATE FOR MISSED SAMPLES 


The concept of sampling assumes the samples are randomly distributed 
such that they faithfully represent the parent population. Some distortion 

wiljl be introduced if the missed samples are not randomly missed. Because 

i ■ ■ . ' 

cloud conditions are the major factors in missed samples, it is unlikely 
that the samples will be randomly missed. It is for this reason that 
this investigation includes regional cloud statistics. Based upon the 
resolution of the simulation, no statistically significant variations with- 
in a region could be detected. This is not surprising for the resolution 
of the cloud model used. Statistical variations from region-to-region were 
evident and were considered when determining the loading effects of over- 
sampling. Before any oversampling procedures could be implemented, it 
is necessary to ascertain the intra- regional biases in the missed samples. 

This could be done with additional simulations on a smal ler scale with 
regionally refined cioud data. The data from this study Indicates a potential 
to reduce the scope of the study to selected regiohs. 






.-Ji 


Another consideration in applying oversampling is the need for multi*’ 
temporal classification. This concept implies that it is not possible to 
cljassify accurately from a single observation. Thus, multiple observations 
wojuld be required for each sany>le and would necessitate an even greater 
Ov|erage for each subsequent observation. The amount of oversampling is a 
function of the percent missed raised to the power of the number of observa- 
tions required. This was not included in the loading consideration of this 
study since only one observation per window was used. 


A third factor to be considered in designing an oversampling strategy is 
the desired confidence level for obtaining samples. For thl^^ study the design 
goal was to achieve the required samples 19 out of 20 years or a 95% confidence 
levei. 


Figure 18 includes the data volumes within the agricultural processing 
system when oversampling is used. The oversampling numbers were determined 
using the regionally weighted data obtained by simulation. The number of 
extra samples required for each region was determined based upon the mean plus 
1.85 sigma of those samples missed in each region. 


SCENE CLOUD EDITING 


The concerns of bias resulting from oversampling lead to a desire to 
minimize the oversampl ing requirement. The oversampling percentages are 
Indicated for each line in Figure 18. The dramatic reduction in the percentages 
Is evident when the samples are extracted from up to 90% cloudy scenes. This 
indicates a very real trade-off in processing costs versus either additional 
satellites or oversampl ing bias. The processing penalty is in the data volume . 
fo be preprocessed. Unfortunately* the limits on the amount of oversampling 
permitted without excessive bias are not currently available to guide the 
assessment of the worth of the additionai processing. Fortunately, as the 
processing Is moved closer to the source, including the ultimate of onboard, 
the vehicles, the concept of editing cloudy scenes vanishes with the entire 
scene concept. 
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COMPARISON OF PROCESSING COSTS 


The previous discussion focused on processing variations in the maximum 
percent of cloud cover acceptable in a scene for preprocessing and the extrac- 
tion of agricultural samples. Some dollar estimates are presented in this 
paragraph as a means for comparing the relative costs of each approach. The 
basis for comparison is the peak processing required during any two week 
period. Peak load impacts the processing system sizing and cost. The peak 
two week load was determined for the two 70k KM satelHte acquisition period 
and the subsequent boards in successive downstream processes were determined 
for three cloud editing strategies. 

The rational used was to develop a system concept for the nominal 
case and normalize the gross order of magnitude cost on a per magapixel basis. 
The acceptance of 50% cloud covered scenes for preprocessing was considered 
nominal . 

The relative processing costs for the 30, 50 and 90 percent cloud cover 
acceptance is illustrated in Figure 19. The agricultural processing is shown 
for the situation of processing all usable information through the models 
and only one sample set per window. In all cases the samples are square 1,2 : ^ 
KM on a side. 

The relative performance of each startegy is measured as a percentage 
of targeted samples obtained. For a performance improvement from 9 A to 99 
percent, the relative processing cost increases from 236K to 517K. This 
provides a strong argument to evaluate the efficacy of oversampling. 

The relative processing costs for each of the editing strategies and over- 
sampling are portrayed in Figure 20. The results of processing oversamples 
^o not appreciably affect the processing cost. By the measure of performance 
tjised, each strategy is equal. The basing effects of oversampling needs to be 
further investigated before a conclusion on the most desirable editing strategy 
can be stated. 
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FIGURE 19. data flow 6 COST COMPARISON 
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AN AGRICULTURAL DATA SYSTEM 


I A simplified concept for an agricultural system employing processing on- 
board a dedicated agricultural satellite is portrayed in Figure 21. This 
system concept is suggested by the results of this study. It is presented as 
a concept for future analysis and has not been subjected to detailed scrutiny. 
Thd availability of the collateral data, the collateral data acquisition 
system, the size and the content of the agricultural database, and the nature 
of the agricultural models and feature extraction algorithms have not been 
unambigously defined. 

The significant findings of this study do Indicate the viability of this 
concept. The nominal Landsat-D orbit appears to be the optiraa'l choice. The 
time of the equatorial nodal crossing needs additional study. 

The significant reduction in data communications and processing require- 
ments by early sample extraction is Justification for considering on-board 
ex^traction. This reduction in the data transmission requirement will make a 
wi'der swath than I 85 KM feasible. The precise orbital and platform stability 
avjai labile with the Landsat-D series of spacecraft will permit registration 
ori-board by means of predetermined mission timelines. The ability to update 
the sample definition is essential consi dering that the sampl ing strategy will 
evolve and can only be optimized by the use of operational system. 

An important simplification of the ground processing will result from 
a { reformatting of the data into sample segments. Subsequent lines of scan 
ddta will contain portions of several segments. The development of sufficient 
o^- board buffers to accumulate an entire sample segment is a practical solution 
to reduce the Identification process. Capacity for several samples to be 
constructed simultaneously will be required. 

On-board radiometric correction is relatively straight-forward. It 
should be done on a scan-by-scan basis prior to the formatting into sample 
segments. It will also be simple to incorporate a cloud editing check since 
tfie acceptance criteria at the sample segment level is either clear or cloudy. 
One-hundred percent cloudy' samples woul d be detectable by the brightness 
threshold..- : 
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AGRICULTURAL SYSTEM 



FIGURE 21. SYSTEM CONCEPT FOR AGRICULTURAL APPLICATION 















It would also be desirable to incorporate some by-pass logic to down- 
link selected portions of the data at various stages in the on-board process 
for engineering verification at the OCC. This system would be used selectively 
for calibrating purposes and would not be expected to impose any additional 
requirements on the communication channel. 

The suggested system concept would use supplemental data from the 
planned Landsat-D system. One use of the supplemental data is for the 
inventory of all land to detect the introduction or loss of agricultural 
land. This could be on an infrequent period of once every 2 to 5 years, 
bamples extracted from the standard Landsat-D processing system could also be 
used to supplement the sample acquisition as well as for episode assessment. 

The size of these samples may be different from those extracted on-board the 
dedicated agricultural satellite. The extraction of agricultural data would 

be done prior to the app 11 cation of any geometric correction resampl ing. 

' ' ' ' , ^ 

k SUMMARY 


When the problems and cost of data processing are factored into the 
results from the data acquisition investigation, one conflict is apparent. 
Better successes are achieved by relaxing the cloud conditions acceptable for 
preprocessing. Lower costs of processing are achieved by making the acceptable 
cloud conditions more stringent. This highl ights the need for more detailed 
stiidies in this area of conf 1 i ct. Some potential resolutions, of which the 
fejjsibiilty is still in question, are the use of oversampling or floating 
sampling techniques and dedicated satel 1 i tes wi th special sensors and on- 
board processing. The costs of development , as well as the techni cal 
practical ityj must be included any future study of these alternatives. 
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APPENDIX A 

ALTERNATIVES TO ADDITIONAL SATELLITES FOR ACQUIRING THE NEEDED INFORMATION 
WITH A MINIMAL LOADING ON THE GROUND PROCESSING SYSTEM 

[ Many factors in a satellite borne image data system affect the number of 
sjatellites required and the amount of ground processing required. The factors 
are coupled such that a change in one parameter will also require a change in 
cither parameters. A qualitative analysis of factors impacting either the 

I 

number of satellites or the amount of ground processing required to perform 
the agricultural mission was performed. The major factors are discussed below. 
-They were considered In baselining the candidate systems investigated in this 
study. 

RELAX REQUIREMENTS 


The requirements for agricultural image data comprise; 1) the crop species, 
2) the makeup of the regions for which the forecast is required, 3) the use of 
the data for inventory, yield, or episode assessment, and k) data quality which 
includes the number of spectral bands, spatial resolutions and spectral resolution 

The crop species of interest, including confuser crops, affect the fre- 
quency and the time criticality of the observations. Relaxation of this 
requirement would have more impact on the volume of data to be processed than 
on the number of satellites required. For the four crops investigated (wheat, 
soybeans, corn, and rice), no impact could be determined on the number of satel- 
lites required. The most significant influence crop species has is its effect 
on the duration of the window during which observations are made. This in turn 
affects discrimination accuracy. This factor is closely related to the data 
quality requirement. The deletion of crops that grow in diverse climates such 
as tropic, artic, or swampland has a slight impact on the total area involved. 

The makeup of the reporting regions has a significant impact on the 
amount of data required to attain a given reporting accuracy. As Castruccio 
expla ins , number of samples required for a region is a function of accuracy 


and not the size of the region. For some regions it is inherently more 
difficult to obtain the required images because of climatological and geo- 
graphical considerations. The orbital parameters are such that regions 
further removed from equatorial latitudes have more frequent opportuni t ies^^^^^ 
observations due to the overlap in the sensor f leld-of-view from orbit-to- 
orbit. Difficult regions were determined using simulation. 

I The intended use of the image data impacts the number of satellites 

I 

required. The use of remotely sensed image data for inventory of the amount 
o|f land planted in a particular crop is easier than to determine the production 
of a crop before harvest. In the former case, only enough observations are 
required to classify the data.* The greatest influence by the use of the 
data, is on the timing of the observation, or window length. The windows used 
for this study were based on a mensuration requirement. In this case, the 
shortest duration of opportunity was 17 days. When yield and episode require- 
ments are added, some of- the windows may be as short as l»-days. While this 
situation was not determined as the baseline for the simulation, some conside- 
rations were given to accommodating such requirements in a maturing system. 

Data quality was not specifically addressed in this study. The general 
relationship of data quality to the data volume and the number of satellites 
required was analyzed. Some generally accepted values were assumed. The 
number of spectral bands affects the data volume in a deterministic way. 

The spatial and spectral resolution requirement affect the number of satellites 
needed in a rather complex relationship of basic sensor response, integration 
tijnes, and the area over which the light is collected. The number of satellites 

'I 

required is dependent upon the ability to observe a particular location at a 
paj-ticular -time. - -As the.swath Js- increased, more locations are in a favorable 
position at a given time. But the^^^^s^^^^ can be increased either by enlarging 
the spatial resolution, bigger steps between integration, or decreasing the 
In teg rat 1 on Vt I me* 


*Thi s statement i s not en t i rel y true as i t app 1 i es to sensor resol ut i on , but 
for this study, existing sensor systems were assumed. 


The latter decreases spectral resolution. For this study, 128 levels of 30 
meter IFQV was assumed. 

ALTER ORBIT 

Alterations in the orbits of the satellites impact the number of satel- 
lites required. The alterations may take the form of; l) frequency of return 
over a particular location, 2) inclination which affects the percentage of time 
over crop regions, and 3) swa^thing which affects the sequence of orbit-to-orbit 
and day-to-day locations of the satel 1 i te. 

The frequency of the return,or repeat cycle, Is a function of altitude and 
inclination. If 100% coverage of the world is assumed, then the number of 
orbits required depends on the swath width of the sensor which was nominally 
assumed at 100 nautical miles. This is the swath width of the currently planned 
Thematic Mapper which Is scheduled for use on Landsat-D. The lower altitudes 
have a shorter period and for the same number of orbits will repeat in fewer 
days. Altitude Is 1 imi ted ip=~the low value by atmospheric drag and the angle 
required wi th the earth surface when the-sensor is at the edge of a swath. 

T^ angle of i 1 lumination of the cropland by the sun during data taking 
affects the classification algorithms. For simplicity a constant sun angle is 
desired and it is obtainable by a proper choice of incl ination 
selected orbital period. Since even for non sun-synchronous inclinations, half 
of each orbital period would not be illuminated and thus only be suitable for 
limited Infrared data, no great advantage was determined for 
orbits.' The detailed investigation of orbital effects using simulation was 
confincled to sun-synchronous orbits. 

Swathing variations include the time sequence of adjacent swaths and any 
uneven distribution of theswaths. The time sequence is particularly interesting; 
as It related to the conditional nature of cloud conditions. When a cloudy 
region is encountered at time T-zeroj there is a high probability a cloudy 
region will be encountered adjacent to that region at T-one if the adjacent 
region is close and if T-pne is not tgo long after T-^zero. For most of the 




ojrbits of interest to agriculture, the time between orbits is of the, order 
ojf one and one-half hours. Likewise, when clear observations are attainable 
ajt one time and place, they are also likely at other nearby places within a 
s^hort period of time. While the net effect may not completely cancel, it is . 
not likely to materially affect the number of satellites required. Of greater 
interest for the agricultural missions is the capability to obtain repeat 
observations of particular areas at selected times. Such effects were studied 
using simulations of shbrt repeat cycle orbits; specifically, 2-day and 9~day. 


USE A PQ INTABLE SENSOR 

j The ability to point the sensor at a specific Sira rget greatly Increases 

* 

tie opportunities to obtain needed information. The amount of pointing consi- 
dered ,ln this study was limited to one swath width to avoid the complications 
0| classifying data obtained at a low grazing angle with the earth's surface. 
Along track pointing also gives an additional dimension for classification 
by taking advantage of anatropic signature discrimination. However, the 
complexity of this relationship to the number of satel 1 i tes requi red was beyond 
t^e scope of this i nvestigatlon. Any system or component requi r ing -redes ign: , 
was rej ected for de ta i led cons i de ra t i on unless the benef its we re outs tand ing. — 
Consequently, the effects of sensor pointing were llml ted to cross-track 
positioning of the currentTy designed Thematic Mapper. 


INCREASE THE SENSOR SWATH 


The swath width used bears directly on the number of satellites required. 
Wilder swaths permit greater spacing between swaths, which results in fewer 
orbits to cover the earth. Thus, more frequent observation opportunities are 
available for a given number of satellites. Any significant change in the 
svath from the 185 KM of the Thematic Mapper will require a sensor redesign. 
Either more data is required or less spatial resolution is obtainable. The 
consideration of swath width is related to spatial and spectral resolution 
discussed under requi rements« Alternatives of increasing the number of sensor 
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primary elements and sensor complexity we're not considered viable. Simulation 
studies of minor swath Increases were performed, A 20% increase In swath width 
(to 222 KM) should be achievable with only minor changes in the present system. 


USE BETTER GROUND TRUTH 


j Better ground truth was considered to have a bearing upon the number of 
satellites required as It relates to the number of observations required to 
obtain sufficient classification accuracy. A consideration was the inclusion 
of an area of high spatlai resolution within each image. The resolving power 
would be Increased optically. An example might be an area 90 meters by 90 
njeters within each 185 KM square Image. The resolution of each picture element 
olf pixel Was assumed at 30 meters. Within the 90 meter squares, one meter 
spatial resolution might be used to establish a training set for classification 
ojf the larger image. The assumption being that the increased resolution 
(pne meter over 30 meter) would permit an accurate classification of the 90 
square meter area. This area would serve as a training set for the nine 30 
meter pixels taken for that area which would be coincident In time and within 
the spatial confines of the image. Insufficient data was available to assess 
quantitatively the effect this would have on the -number of-satel lltes required.. 
Once the effect on classification Is determined j the results obtained for the 
assumed classlflcat^^ requirements can be sea led. | 


USE SAMPLING STRATEGIES 


j The effective use of sampling strategies achieves both a reduction In the 
anrount of data to be processed and the number of satellites required. Only a 
Small percentage of the data obtained will be used In the agricultural fore- 
casting models. The selection of the right data by sampling at the data taking 
phase effectively reduces the data. The statistical nature of the data and 
the chance occurrence that the data will not be rendered useless because of 
cloud cover places an increased demand on the number of observation oppor 
tuni ties required. An acceptance of the condition that less than all the 
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desired samples will be obtained reduces the pressure for more satellites. 

Then, an oversampling strategy can be employed such that an attempt is made 
to obtain more samples than are actually needed. This requires the assump- 
tion that the designated samples contain equally valid information. This 
1 limits the application to some small percentage of the total samples. For 
an arbitrary limit of success can be claimed in meeting the requirement 

wijenever the number of satellites employed obtained 90i of the targets specified. 

,, Tnus, if SB% of the samples were required to achieve a given mensuration 
accuracy, and the simulation showed 92% obtainable with two satellites, over- 
simpling strategies would call for an attempt at acquiring 108.7% of the 
number of desired samples. The penalty Is an increase in the amount of data 
to be processed for the benefit of requiring fewer satellites. 

PRIORITIZE DATA ACQUISITION 

I The nature of the agricultural mission places a temporal value on the 
i|)formation obtained. There is also some value to second and third observa- 
tions of basically the same images because they help reduce classification 
errors. There is a trade-off between redundant observations and the cost of 
additional capabl 1 ity to process them. A priority strategy will maximize 
the information throughput for a given processing capacity. ^ ^ 

' - The use of a priority system for data acquisition indicates a 
coupling between the acquisition subsystem and the user models since the data 
value will dynamically change according to user model needs, 

EDIT ON IMAGES 

Certain preprocessing, which as a minimum, includes the correlation of 
the image data to time and position via spacecraft engineering data, must be 
performed on all usable data. Radiometric correction may also be included 
in the preprocessing. To minimize the amount of unproductive preprocessing, 
some early Image editing Is desirable. The exclusion of 100% cloud covered 
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images Is a reasonable way to reduce the amount of processing performed on 

data with little likelihood of producing usable Information. As the quality 

1 

requirements of a scene acceptable for processing Is Increased, the propor- 
tion of useless processing decreases. The overall success In acquiring a 
specified number of samples for a fixed number of opportunities, (as a function 
of the number of satellites used) will also diminish. The quantifications of 
these effects was done using simulation. 


INHIBIT ACQUISITION OF POOR QUALITY DATA 


, The development of a sensor system to perform editing at the point of origin 
wl p reduce the processing even more. The approach to this alternative Is to 
quantify the reduction obtainable at each point In the system wl thout regard 
to the technology or cost of performing the reduction. The resulting benefits 
will provide a potential measure of the worth of alternative systems. 

USE STABLE PLATFORM 

The use of an ultra stable platform to minimize regl strat Ion -dl f f Icul t les and 
the amount of geometric correction was considered. The current projection for 
Lahdsat-D when combined with about a one KM square sample segment were considered 
justification for not considering geometric correction as a necessary part of 
agricultural processlng.^^^^^^ 


IMPROVE PROCESSING 

An improvement In processing Is a definite alternative to the problem of 
editing and ground processing. It also can take various forms ranging from a 
reduction in the number and frequency of inputs required to the models to 
more efficient classification algorithms. While admitted as an alternative with 
great potential. It was deleted from consideration In this study. Any Improve- 
ment will benef I t the system in addition to the recommendations from this study, 
without rendering the study eccentric to the other alternatives considered. 
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APPENDIX B 


MAXIMUM SWATH WIDTH AS A FUNCTION OF ALTITUDE 

Precise information on the maximum off~nadlr angle that is acceptable for agri- 
cultural images is not avai labie. Some gui del ines are extrapolated from experience 
with Landsats i and 2. For a nominal 900 KM altitude the maximum usable off-nadir 
angle was found to be 20 degrees. 


Earth Radius = 6378. 16$ KM 



Solving the two equations for X and y yields 
X = 330.6956 

y - 6369.58251 

-'"'V ■ ■■■ ■ J 

This yields a solution of 6 m tan 2.97201° 


The mi nimijm acceptable grazing anjgle3 is then found by the geometry of the sum 
of the interior angles of a triangle equaTJ 

$ = 6 * -90 where a + 6 +3 ' = 180 

or 20 + 2.97201 = 180 

3 = 180° - 20° - 2.97201° -90° =67.03 

To obtain the general relationship of altitude as a function of 6 as constrained 
by a minimum 3 of 67.03°» it Is necessary to apply the sine law: 

• A b: £ • ■■ ■■ 

sin ot " sing SinJ^ 

Then radius of earth + altitude = line A-B — radius of earth 
sin (90+67.03°) Sin 0 Sin a 

App lying the constants yeilds 

6378.165 + h = S = 6378.165 
0.3902^*90997 Sin 6 Sin a 

Apply the constraint that 3^ + 0+ a= 180° which yields a* (22.97 ■ 0 ) 


The exp ress ion then becomes 


6378.165 + h = 6378.165 ... . 

0.39Q2490997 Sin (22.97- 6) reduces to 


h = 21*89.0731^9 
Sin (22.97 -0) 


- 6378.165 
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APPENDIX C 


INVESTIGATION PROCESS 

The sequehce of steps involved in investigating each of the candidate 
systems is illustrated in Figure C-1. 

STEP 1 - For each of the candidate systems, the orbit parameters were 
calculated using Reference 5- 

STEP 2 - It was necessary to alter the nominal orbit parameters 

[altitude and inclination) slightly to obtain a repetition of 
i:he nadir equational crossing after the proper number of 
orbits. This was necessary because a one degree variation 
o£ the end of the repeat cycle is equivalent to a 111.32 KM 
error in the repeat swath at the equator. For all orbits 

I 

exciept the ll»85 KM orbit, the swath ground trace repeated with 
.7^ KM. For economy of computer run time, the DSDS Mission 
Ephemeris Generator (MEG) was only run for one repeat cycle 
for each latitude and insertion point. The mission ephemerides, 
in terms of nadir latitude and iongitude, were generated in along 
track increments of 90 nautical miles, which corresponds to 
the along track spacing of standard Landsat scenes. This 
permitted a direct accounting of scenes as a measure of data 
. ' volume suitable for comparison with previous studies. 

STEP 3 - The data generated by MEG was used to generate pairs of 

latitude and longitude corresponding to the edge of the sensor 
swath. For this study , nadi r pointing was usei exGTusively. 

For increased swath width^^^^^ it was not necessary 

to repeat the MEG runs. 

STEP 4 The Target Model was run for each candidate system comprising 
a combination of satellites at different altitudes, insertion 
points and swath widths. This model determines, by target 
number, which targets are within the sensor swath for each 
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FIGURE C-l. INVESTIGATION PROCESS 










step (scene) during the complete repeat cycle. This data 
is saved for repeated use as input to the crop model. A typical 
report from the Target Model is shown in Figure C-2. The report 
shows how many times each target was seen during a repeat cycle 
and summarizes the results. As seen in Figure C-2, some targets 
are observed several times during a repeat cycle, This is due 
to overlap between adjacent swaths which, for a I 85 KM swath, 
reaches 100^ at 57.^6° latitude. 

STEP 5 “ the final step in the investigation of the candidate 

orbits involved running the Crop Model. For input data, the Crop 
Model will accept up to three Target Model output files, 
the following information required by the Crop Model is input 
on parameter cards: 

o Simulation Start Day 
o Simulation End Day 

0 Number of observations of a target to be processed during 
(A record is kept of the number of viewing 
opportunities and the number of opportunities that are 
^ free. Only the number specified are placed in the 

time line file for processing) 
o Level of scene cloud cover acceptable for processing 
o the number of satellites, and for each satellite; 

- Al ti tude 

- Inclination 

- Number of orbits per repeat cycle 

- Period per orbit in seconds 

[■ For each test case, the simulation was a ful 1 year. When the end 
of the repeat cycle was reached on any of the input files, the file was re- 
wound and reread until the full year was covered. 

After a scene record was read from each satel 1 i te ' s input file, the 
following steps were performed. 

A. A test was made to determine which scene occurred first. 

B. A' check was made to determine if any of the targets in the 

' scene had crops in a growth stage of interest , i .e. , 

an active window. If none of the targets were active, the next 
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TARGET NUMBER OF OBSERVATIONS 
INDEX NUMBER DURING REPEAT CYCLE 
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FIGURE C-2. typical OUTPUT FROM THE TARGET MODEL 


scene was read from the corresponding input file and Step A was 
repeated. 

C. If any of the targets in the scene were in an active window, the | 
Cloud Model was called to determine the cloud category for the scene 
based on the current month, time of day and the cloud region. 

D. For each target in an active window, the number of opportunities 
for observation was incremented. 

E. Thje scene cloud cover was compared with the scene cloud cover 

threshold for processing. If the scene cloud cover was below 
thL threshold, a test was made on each target in the scene. The 
percent of targets that could be obtained from a scene is a function 
of the scene cloud cover as shown below. ^ ' 


Scene Cloud 
Category 


Percent of 
Cloud Cover 
in the Scene 


Percent of Clear 
Samples lObtai'nable 


1 V 10 • 95- 

2 10, 20, 30 80 

3 M, 50 55 

k 60, 70, 80, 90 25 

5 100 0 

A separate test was performed for each target. if the target was^ 
clear, the number of clear observations was incremented. j 

F. The Target Processing Timeline was created. For each clear observation 
a] test was made to see if the target had been, seen for the desired 
niiimber of times during the current window, if the desired number 
of clear observations had not been obtained, a record was written 
in the t ime 1 ine file with -the following information: 


Target number 
Date 

Crdp(s) t>f interest 
Window number 



^ Scene Processing Timel ine was created. For each Scene with at 

least one target in an active window, the following information was 
written into the time line file: 

Date 

- Region number 
Cloud region 

' ~ Scene cloud category 

At the end of the run, the Crop Model report generator was called. 
Typical outputs from the Crop Model report generators are described 

;;y; ■ , ":next.r " 

CROP MODEL REPORTS 

The next four figures show a portion of the reports’ from a typical 
Xrop Model run. The first listing, shown in Figure C~3, contains the input 
parameters for the run. The case shown is for a one year run with 2 Landsat^D 
satellites;' In this run, all scenes with 50% ploud cover or less were accepted 
for processing. The maximum number oJ' observations to be processed for 
anf Carget during a window was set at 2. - 

The scenes and sample segments acquired for processing on a daily basis 
are presented in Trade Study Report 1, A portion of Report 1 is shown in 
Figure C-A. The day of the year and date are given in the first 3 columns. 

The number of point targets (sample segments) to be processed for the day is 
listed in Column A. For an operational system, this number would be scaled 
up by a factor of approximately 30, The next 5 columns show the number of 
scenes acquired in each cloud category. A scene is recorded only if there is 
« at least one target in an active window within the ssnSor swath. The last 

three columns show the total scenes acquired for "the day)), the number of scenes 
with acceptable cloud cover and the percent of the total^*lcenes t^ were 
acceptable. The last 1 ine of Trade Study Report 1 gjyes the yearly totals. 

For the week starting August 27, as shown in Figure C-l», there were a total 
of 517 scenes acquired of which 337 had 50% cloud cover or less. From these 
scenes, 159 sample segments were extracted for processing. For an operational 
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THE INPUT PARAMETERS hOK THIS RUN ARE J 

1 = SIMULATION START DAY 
'365 = SIMULATION END DAY 

233 = NUM3ER OF ORBITS IN A REPEAT CYCLE FOR VEHICLE NUMBER I 
5933 = NUMBER OF SECONDS IN AN ORBIT FOR VEHICLE NUMBER I 
98.20 = ORBIT INCLINATION IN DEGREES FOR VEHICLE NUMBER 1 
704 = ORBIT ALTITUDE IN KM. FOR VEHICLE NUMBER 1 
185 = SWATH WIDTH IN KM. FOR VEHICLE NUMBER I 
233 = NUMBER OF ORBITS IN A REPEAT CYCLE FOR VEHICLE NUMBER 2 
5933 = NUMBER OF SECONDS IN AN ORBIT FOR VEHICLE NUMBER 2 
98.20 - ORBIT INCLINATION IN DEGREES FOR VEHICLE NUMBER 2 
704 = ORBIT ALTITUDE IN KM. FUR VEHICLE NUMBER 2 
185 = SWATH WIDTH IN KM. FOR VEHICLE NUMBER 2 

3 * MAXIMUM CLOUD CATEGORY ACCEPTED F OR PROCESSING 

2 * MAXIMUM NUMBER OF OBS. PROCESSED DURING A WINDOW 


Ft GURE C- 3 . L I ST I NG OF I NPUT PARAMETERS FOR A CROP MODEL RUN 


^^Rqvauty 
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FIGURE C-4. SAMPLE SEGMENT AND SCENE ACQUISITION TIME LINE 



system, this would equate to |59 x 30 = i»770 sample segments and represents 
a processing load of 95^ sample segments per day for a 5"day week. 

A summary of the scenes observed for the year In each region is shown in 
Figure C- 5. The information presented is the same as that contained in the 
last -8 columns of Trade Study Report 1. There were no scenes acquired for 
Region 31 (USA E) because none of the four major crops grew in sufficient j 

quantities In this region; therefore, no targets were assigned to the region. 

The percent of scenes with 50% cloud cover or less ranged from a low of kk.7% 
in the central region of India to a high of 84.8% in Egypt. 

1 ','7' A 

^ C-6 is typical of the reports for each phenol og I cal region. For each 

tar 9 ®t assigned to the region, the following information is given: 

o f|heno1ogical Region Number - The 36jg®pgraphic r^ on Table 2, page i 

divided Into 58 homogeneous pheno log leal regions, and 
wi th a unique crop calendar. 

0 Target Number - Target numbers were assigned consecutively In each 
country. 

o Cloud Region - Specif ies-which of the 30 cloud model regions the 
' target is in. ■ ^ , v 

o latitude and Longitude. 

^ ^ during which observations are to be obtained, five 

additional pieces of information are given: 

Crop Code A one in any of the five columns under the crop code 
means that an image during the window is used to discriminate that 
pi^op from other crops in the region. (The crop code is defined 
^-r ^ second line of the report.) 

- Start - Is the day of the year on which the window starts. 

- LEN - Is the number of days that the window is active beyond the 
start -day. 

- 08s ’• Is the number of times that the target was within the swath 
of the satel lite's sensor during the window, i.e. , the number Of 
observation opportunities regardless of cloud cover. 

*These phenological regions are regions with equivalent climates, soil con- 
ditions, and cropping practices that may be characterized by a single crop 
calendar. • 



REGION 


SCENES OBSERVED BY 

CLOUD CATEGORY 

TOTAL 

SCENES 

ACCEPTED 

: PERCENT 
ACC FOR 

NUMBER 

1 NAME 

CC 

1 CC 2 

0*^0 '3 

CC 4 

CC 5 

SCENES 

FOR PROC 

PROCESSING 

1 r 

ARGENTNA 

7^3 

105 

33 

85 

17 

313 

211 

67.4 “ 

- 2 

AUStRLIA 

50 

49 

15 

40 

6 

160 

il4 

71.2 

3 

0 ANGLOS H 

9 

24 

15 

26 

5 

79 

48 

60.8 

4 

BRAZIL N 

36 

62 

40 

79 

9 

226 

138 

61.1 

5 

BRAZIL S 

46 

76 " 

35 

103 

23 

283 

157 

55.5 

6 

CANADA 

38 

141 

39 

126 

18 

362 

218 

60.2 

7 

CHINA N 

37 

169 

45 

204 

28 

483 

251 

52.0 

. 8 

CHINA C 

i53 

l79 

55 

170 

39 

604 

387 

64.1 

9 

CHINA S 

51 

66 

21 

65 

12 

215 

138 

64.2 

10 

1 EGYPT 

42 

35 

18 

15 

2 

112 

95 

84.8 

11 

FRANCE 

60 

85 

75 

95 

24 

339 

220 

64.9 ^ 

12 : 

INO PUNJ 

28 

37 

15 

21 

3 

104 

80 

76.9 

13 

In6 GANG 

lu 

29 

19 

37 

12 

107 

58 

S'*. 2 -0 a 

^‘‘•2 §3 

54.0 wg 

54.1 of; 

14 

INO CENT 

20 

37 

31 

85 

24 

197 

88 

15 

INO BILA. 

15 

25 

21 

40 

12 

113 

61 

16 

INO COST 

30 

53 

36 

82 

19 

220 

119 

17 

INOONESA 

30 

54 

33 

49 

7 

173 

117 

67.6 

13 

ITALY 

70 

75 

69 

74 * 

11 

299 

214 

71.6 ^ Q 

19 

JAPAN 

18 

25 

27 

38 

7 

115 

70 ■ 

■■”60.9-— 3^ 

20 

MEXICO 

92 

160 

60 

162 

40 

534 

332 

62.2 ® 

21 

PAKISTAN 

23 

32 

8 

24 

4 

91 

63 

69.2 

22 

RUMANIA 

41 

47 

27 

39 

11 

165 

115 

69.7 

23 

S.AFR-ICA 

26 

!! 35^ ■ 

18 

40 

13 

132 

79 

■ 59.8 

24 

PHILPNES 

44i 

•i ■■ 

26 

74 

18 

203 

111 

54.7 

25 

THAKLANO 

13 

36 

31 

60 

1 1- 

151 

80 ■ 

■ 53.0 

.26 _ 

TURKEY 

31 

22 

45 

28 

3 

129 

98 

76.0 

27 

'‘ U.S. A 

69 

84 

49 

97 

21 

320 

202 " 

63^1 

28 

U.S. B 

78 

93 

80 

101 

3Z 

384 

251 

65.4 

29 

U.S. C 

130 

164 

125 ^ 

159 

35 

613 

419 ■ ■■*■* 

68.4 

30 

U.S. 0 

94 

119 

42 

142 

35 

432 

255 

59.0 

31 

U.S. E 

0 

0 

0 

0 

0 

0 

0 

C.O 

32 

USSR LAT 

134 

^ 167 

127 

152 

38 

618 

428 

69.3 

33 

USSR UKR 

25 

51 

29 

48 

10" 

' 163 

105 

64.4 

34 

USSR T-V 

116 

106 

101 

95 

17 

435 

323 

74.3 

35 

USSR V-U 

92 

153 

94 

168 

39 

546 

339 

■ “62.1* ■ 

36 

USSR SIB 

47 

158 

47 

139 

23 

414 

252 

60.9 

37 

YUGOSLAV 

41 

58 

41 

45 

17 

202 

140 

69.3 


FIGURE C-5. 

SUMMARl 

r OR SCENES 

OBSERVED 

BY REGION 

AND CLOUD 

CATEGORY 

FOR THE YEAR 





TRADE STUDY REPORT 2 -POINT TARGET DAT A 
FRANCE FRANCE FRANCE FRANCE CROP CuOE = iWCSR /sSpIung WHEAT RfWiNTER WHEAT C»C0RN,. S»S0VBEANS^^R«RICE. 









WI NUUW 

1 




WINDOW 2. _ 


WINDOW 3 

PHENd 

TARGET 

CLOUD 



CHOP 





percents 

CROk 



PERCENT 

CROP PERCENT 

REGION 

NUMBER 

REGION 

LATITUDE LONG 

CODE 

STAkT 

LEN 

OrtS CIR 

CLEAR 

CODE 

START 

LEN OBS^CLR 

JCLEAR _ 

_.CO3E_STA.RJ_LEN_O0.S_CJJL-CLEAR 




: - 


ZMCSR 






ZWCSR 




ZWCSR 

~ 11 ” 

4Z9 

ll 

13 a . 08 

4.30 

1000 

120 

31 

6 

2 

33.3 

100 

162 

59 12“ 4” 

33.3 - 


. 11 

V30 

II: 

133.14 

359.97 

1000 

1 20 

31 

4 

1 

25.0 

100 

162 

59 8 4 

50.0 


11 

All 

^ ■' ll: ■■ 

133.63 

4. 90 

1000 

120 

31 

6 

5 

83.3 

100 

162 

59 12 7 

58.3 


. ll 

«2 

M 

139.81 

2.47 

1000 

120 

31 

6 

2 

33.3 

ICO 

162 

59 12 9 

75.0 


11 

411 

11 

137.70 

358.17 

1000 

120 

31 

6 

1 

10.7 

100 

162 

59 10 5^ 

50.0 


. . 11 

_ 414 

11 

140.80 

2.79 

1000 

120 

31 

6 

3 

5C.0 

100 

162 

59 12 6 

50.0 


11 

435 

11 

135.27 

2.88 

1000 

1 20 

31 

a 

4 

50.0 

100 

162 

59 16 9' 

56.3 


_ ll 

_ 436 

11 

133.84 

2.92 

1000 

120 

31 

4 

1 

25.0 

100 

162 

59 8 3 

37.5 


11 

• 43T 

ll 

139.88 

4.14 

1000 

120 

31 

6 

4 

66.7 

100 

162 

59 12 4 

33.3 


^ ll _ 

^ A18 

11 

135.62 

359.63 

1000 

120 

31 

3 

2 

2 5.0 

100 

162 

59 16 8 

50.0 


11 

439 

11 

139.77 

1.60 

1000 

120 

31 

a 

. 5 

62.5 

100 

162 

59 16 8 

50.0 


_ 11 _ 

^.440 

ll 

140.14 

2.65 

1000 

120 

31 

6 


66.7 

100 

162 

59 12 7 

58.3 


11 

441 

11 

136.33 

3.07 

1000 

120 

31 

8 

2 

2 5.0 

100 

162 

59 16 8‘ 

.50.0' 


_ 11 

^.442 

11 

138.58 

3.90 

1000 

120 

31 

8 

5 

62.5 

100 

162 

59 16 9 

56.3 


11 

443 

11 

140.69 

1.79 

1000 

120 

31 

8 

3 

3 7.5 

100 

162 

59 16'^ 13“ 

81.3 


.11 ^ 

^44 

iliir. 

134.40 

6.25 

1000 

120 

31 

4 

4 

100.0 

100 

162 

59 8 1 

12.5 


11 / 

445 

ll 

136.48 

0.70 

1000 

120 

31 

4 

1 

25.0 

lod 

162 

59 '^’ 8 4 

50.0 



11 

-_.A46 

H 

134.24 

3.14 

1000 

120 

31 

4 

2 

50.0 

100 

162 

59 8 3 

37.5 


11 

4^7 

11 

138.64 

7.32 

1000’ 

120 

31 

8 

6 

75.0 

' ’ 100 

162 

59 16 id 

62.5 


ll 

_|440 

ll 

139.28 

5.38 

1000 

120 

31 

8 

4 

50.0 

too 

162 

59 16 8 

50.0 

V 

ll 

449 

11 

139.70 

1.83 

1000 

120 

31 

8 

4 

50.9 

100 

162 

59 16 To 

62.5 

-jP : ^ 

^ 11 w 

450 

11 

137.91 

0.32 

1000 

120 

31 

6 

2 

33.3 

100 

162 

59 12 7 

58.3 


11 

451 

, ■‘^ ll ■ 

136.2? 

359.19 

1000 

120 

31 

6 

3 

50.0 

100 

162 

59 11 “7“ 

63.6 

-— — — g 

Hi 

452 

_ 11 

,137.50 

1.08 

1000 

120 

31 

4 

3 

75.0 

LOO 

162 

59 8 3 

37.5 


Ilf 

451 

11 

1140. 31 

3.69 

1000 

120 

31 

8 

3 

37.5 

100 

162 

59 16 6 

37.5 


11 

454 


133.38 

357.22 

1000 

120 

31 

6 

1 

16.7 

100 

162 

59 ll 4 

36.4 


ll 

455 

ll 

139.73 

2.34 

1000 

120 

31 

6 

3 

50.0 

100 

162 

59 12 6^ 

50.0 


__ 11 

456 

11 

133.28 

2.40 

ICOO 

120 

31 

6 

3 

50.0 

100 

162 

59 12 5 

41.7 


11 

457 

ll 

138.45 

3.63 

1000 

120 

31 

a 

2 

25.0 

100 

162 

59 ' 16 8* 

50.0 


^11 

458 

ll 

134.22 

3.91 

1000 

1 20 

31 

8 

2 

25.0 

100 

162 

59 16 10 

62.5 


11 

459 

11 . 

133.49 

'358.87” 

1000 

120 

31 

4 

. 1 ' 

25.0 ‘ 

'lOO 

162 

59' 8 3 

37.5 


11 

^460 ^ 

11 

138.80 

6.91 

1000 

120 

31 

8 

5 

62.5 

100 

162 

59 16 10 

62.5 


ll 

: 461 

11 

134.14 

7.42 

'1000 

120 

31 

6 

3 

50.0 

100 

162 

59 10 "5 

50.0 


^ 11 

462 

ll 

133.24 

0.32 

1000 

120 

31 

6 

2 

33.3 

100 

162 

59 12 6 

50.0 


11 

463 

ll 

133.78 

4.64 

1000 

120 

31 

4 

4 

100.0 

100 

162 

59 0 '* 5 

62.5 


^ll-_ 

_ 464 


13 7.39 

359.01 

1000 

120 

31 

6 

0 

0.0 

100 

162 

59 10 7 

70.0 


ll 

465 

11 

134.96 

0.75 

1000 

120 

31 ' 

6 

4 

66.7 

100 

162 ' 

59 12 6 

50.0’ 


. 11 


11 

135.00 

5.67 

1000 

120 

31 

6 

5 

6'2.5 

100 

162 

59 16 7 

43.8 


11 

467 

M;,. ll: ^ 

139.45 

5.10 

1000 

120 

31 

6 

2 

33.3 

100 

162 

59 12 6 ‘ 

50.0 ' 


u _ 

ll 

139.02 

0.73 

1000 

120 

31 

6 

3 

50.0 

100 

162 

59 12 10 

83.3 


ll 

469 

ll 

139.79 

4.43 

1000 

120 

31 

a 

5 

62.5 

100 

162 

59 16 12" 

75.0 ■' 


ll „ 

470 

11 

137.11 

358.59 

1000 

120 

31 

6 

3 

50.0 

100 

162 

59 10 7 

70.0 


. 11 

471 

11 

136.24 

338.96 

LOGO 

120 

31 

6 

3 

50.0 

100 

162 

59 ll 6 

54.5 


^ ll 

472 

11 

134.99 

3.81 

1000 

120 

31 

6 

3 

50.0 

locr 

162 

59 12 6 

50.0 


11 

473 

ll. 

■\ ' ' V 

134.27 

5.29 

1000 

120 

31 

8 

5 

62.5 

100 

162 

59 16 8 

50.0 ‘ 



: ■ ■ ■ V 

■] ■. 1 : : 

TOTALS 

FOR WINDOW 1;- 
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A5 SAMPLES IN THIS REGION 

0 WINDOWS HAD ZERO OPPORTONi Tits " 

1 WINDOWS HAD ZERO CLOUD iFREE OUSERVATIONS - 

F I GURE C-6 . TYP I CAL REPORT FOR TARGETS I N A PHENOLOG I CAt REGTSN ^ 


- CLR - Is the number of times that a clear (cloud-free) observation 
was obtained during the window. 

- Percent Clear -Is the percent of the total observational 
ojjportuni ties which were clear. 

At the bottom of the report, the number of observational opportunities 
and the number of clear observations are totaled for each window. The percent 
of thd total opportunities that were clear is also given. 

A summary of the number of targets in the region, the number of windows 
that did not achieve at least one observational opportunity and the number of 
windows that did not achieve at least one cloud free observation are also given. 
For the example shown, all targets had observational opportunities and only 
Target Number 47^ failed to have at least one cloud free observation. 


During the first window, shown in Figure C-6, each target was observed 
either k, 6, or 8 times. In this run, there were two satellites, each with a 
l6-day repeat cycle. Since the window has a 32-day duration, each satellite 
wouild see each target twice exclusive of overlap coverage. If a target was in 
the overlap coverage area for one satellite, it would be observed six times, and 
if a target was in the overlap coverage area for both satellites, it would be 
observed eight times. 


STANDARD TEST CONDITIONS 


For all test cases to be discussed, the following conditions are standard 
unless otherwise specified: 


-Swath Width - A standard Landsat sensor with a Swath width of 185 KM i 
(1|00 Nautical Miles) was used. 

- Scene - A scene represents an area one swath width wide with an along 
track distance of 90 Nautical Miles. 

Cloud Cover Accepted for Processing - All scenes with less than or 
equal to 50% cloud cover were accepted by the preprocessor for 
sample segment extraction. 


- Ol^servations Processed per Window - Only one clear observation was 
requi^ during each window. For regions with multiple crops, a 
Separate w defined to discriminate each crop from the other 

confusion crops. 

The various test cases are compared on their capability of obtaining 
needed sample segments* The goal for acquisition of samples was set at 
98% of the designate samples for each region. This level was set to minimize 
the mensuration error and the introduction of bias caused by obtaining a 
disporportionate share of samples from areas that are relatively cloud free. 


APPENDIX D 

CLOUD MODEL 6 STATISTICS 


The Cloud Model used in this simulation is identical to the cioud model 
usied in a previous study, it is described in detail » "Global 
Crop Production Forecasting Trade Study - Volume II - Approach and Results," 
Section If- 14. 


j Tables 1 and 2 present cloud statistics based on the twenty year runs 
with two Landsat-D satellites. Table 1 gives the mean and standard deviation 

I 

fojr the percent of clouds in cloud categories I through 3^ (less than 50 % 
cl'oud cover) for each region. Table 2 contains the mean and standard devia- 
tion for the 20 year period regardless of region. On the average, 63.2% 
of the scenes had 50% cloud cover or less. 
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TABLE 1. MEAN AND STANDARD DEVIATION OF THE 
PERCENT OF SCENES WITH 50!^ CLOUD 
COVER OR LESS 


REGION 

NUMBER 

NAME 

PERCENT OF SCENES 
IN CLOUD CATEGORIES 1-3 

MEAN 

STANDARD DEVIATION 

1 

Argent i na 

67.5 

2.3 

2 

Austral ia 

76.6 

3.2 

3 

Bangladesh 

52.8 

5.5 

k 

Brazil North 

61.7 

3.3 

5 

Brazil South 

57.2 

3.0 

6 

Canada 

55.7 

2.5 

7 

China North 

50.7 

2.5 

8 

China Centrai 

64.0 

2.2 

9 

China South 

66.3 

3.0 

10 

Egypt 

86.8 

3.1 

11 

France 

66.9 

2.8 

12 

India Punjab 

71.4 

4.5 

13 

India Gange ; 

47.8 

6.6 

14 

India Central 

46.2 

3.2 

15 

India Bilaspur 

50.5 

4.6 

16 

India Costal 

58.1 

3.5 

17 

Indonesia 

63.4 

4.2 

18 

Italy 

68.8 

2.6 

19 

Japan 

57.2 

4.9 

20 

Mexico 

64.5 

1.9 

21 

Pakistan 

69.6 

5.3 

22 

Roman i a 

66.8 

3.6 

23 

South Africa 

66.8 

3.6 

24 

Phi 1 ippines 

52.5 

3.1 

25 

Thai land 

57.7 

3.9 

26 

Turkey 

69.1 

3.9 

27 

USA - Region A 

65.7 

3.1 

28 

USA - Region B 

68.1 

1.9 

29 

USA - Region C 

66.6 

2.0 

30 

USA - Region D 

59.5 

2.7 

31 

USA - Region E 



32 

USSR Latvia 

67.7 

2.0 

33 

USSR Ukraine 

67.0 

4.6 

34 

USSR Transvolga 

69.9 

2.3 

35 

USSR Voiga-Ural 

65.9 

1 .4 

36 

USSR Siberia 

56.4 

2.4 

37 

i 

Yugoslavia 

66.3 

2.2 


7 ^» 
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TABLE 2. MEAN AND STANDARD DEVIATION FOR THE 
OCCURRENCE OF EACH CLOUD CATEGORY. 


Cloud j,,'. 

Categoryf i 

Percent of 
Cloud Cover 

PERCENT OF SCENES IN CATEGORY 

Mean 

Standard Deviation 

: :i ■■ 

' 0 

19.5 i 

1 ' 
; 


0.26 

x2\ ' 

10, 20, 30 

27.9 

' j 
• r 

. 

0.35 / 


AO, 50 

15.8 ' 


■ 

: 0.30 


60, 70i 80, 90^ 

30.0 



0.40 

5 : ' ^ 

100 . 

6.8 



0.29 

1 - 3 : 

0 - 50 

63.2 



0.40 






